An optimized approach for multiplexing single-nuclear ATAC-seq using oligonucleotide-conjugated antibodies.

BACKGROUND: Single-cell technologies to analyze transcription and chromatin structure have been widely used in many research areas to reveal the functions and molecular properties of cells at single-cell resolution. Sample multiplexing techniques are valuable when performing single-cell analysis, reducing technical variation and permitting cost efficiencies. Several commercially available methods have been used in many scRNA-seq studies. On the other hand, while several methods have been published, multiplexing techniques for single nuclear assay for transposase-accessible chromatin (snATAC)-seq assays remain under development. We developed a simple nucleus hashing method using oligonucleotide-conjugated antibodies recognizing nuclear pore complex proteins, NuHash, to perform snATAC-seq library preparations by multiplexing. RESULTS: We performed multiplexing snATAC-seq analyses on a mixture of human and mouse cell samples (two samples, 2-plex, and four samples, 4-plex) using NuHash. The analyses on nuclei with at least 10,000 read counts showed that the demultiplexing accuracy of NuHash was high, and only ten out of 9144 nuclei (2-plex) and 150 of 12,208 nuclei (4-plex) had discordant classifications between NuHash demultiplexing and discrimination using reference genome alignments. The differential open chromatin region (OCR) analysis between female and male samples revealed that male-specific OCRs were enriched in chromosome Y (four out of nine). We also found that five female-specific OCRs (20 OCRs) were on chromosome X. A comparative analysis between snATAC-seq and deeply sequenced bulk ATAC-seq on the same samples revealed that the bulk ATAC-seq signal intensity was positively correlated with the number of cell clusters detected in snATAC-seq. Moreover, when we categorized snATAC-seq peaks based on the number of cell clusters in which the peak was present, we observed different distributions over different genomic features between the groups. This result suggests that the peak intensities of bulk ATAC-seq can be used to identify different types of functional loci. CONCLUSIONS: Our multiplexing method using oligo-conjugated anti-nuclear pore complex proteins, NuHash, permits high-accuracy demultiplexing of samples. The NuHash protocol is straightforward, works on frozen samples, and requires no modifications for snATAC-seq library preparation.

Genomic insights into host and parasite interactions during intracellular infection by Toxoplasma gondii.

To gain insights into the molecular interactions of an intracellular pathogen and its host cell, we studied the gene expression and chromatin states of human fibroblasts infected with the Apicomplexan parasite Toxoplasma gondii. We show a striking activation of host cell genes that regulate a number of cellular processes, some of which are protective of the host cell, others likely to be advantageous to the pathogen. The simultaneous capture of host and parasite genomic information allowed us to gain insights into the regulation of the T. gondii genome. We show how chromatin accessibility and transcriptional profiling together permit novel annotation of the parasite's genome, including more accurate mapping of known genes and the identification of new genes and cis-regulatory elements. Motif analysis reveals not only the known T. gondii AP2 transcription factor-binding site but also a previously-undiscovered candidate TATA box-containing motif at one-quarter of promoters. By inferring the transcription factor and upstream cell signaling responses involved in the host cell, we can use genomic information to gain insights into T. gondii's perturbation of host cell physiology. Our resulting model builds on previously-described human host cell signalling responses to T. gondii infection, linked to induction of specific transcription factors, some of which appear to be solely protective of the host cell, others of which appear to be co-opted by the pathogen to enhance its own survival.

ETS1 is a novel transcriptional regulator of adult T-cell leukemia/lymphoma of North American descent.

Adult T-cell leukemia/lymphoma (ATLL) is an aggressive T-cell lymphoma associated with the human T-cell lymphotropic virus type 1 virus endemic in regions including Japan, the Caribbean islands, and Latin America. Although progress has been made to understand the disease, survival outcomes with current standard therapy remain extremely poor particularly in acute ATLL, underlying the need for better understanding of its biology and identification of novel therapeutic targets. Recently, it was demonstrated that ATLL of North American-descendent patients (NA-ATLL) is both clinically and molecularly distinct from Japanese-descendent (J-ATLL), with inferior prognosis and higher incidence of epigenetic-targeting mutations compared with J-ATLL. In this study, combined chromatin accessibility and transcriptomic profiling were used to further understand the key transcriptional regulators of NA-ATLL compared with J-ATLL. The ETS1 motif was found to be enriched in chromatin regions that were differentially open in NA-ATLL, whereas the AP1/IRF4 motifs were enriched in chromatin regions more open in J-ATLL. ETS1 expression was markedly elevated in NA-ATLL in both cell line and primary tumor samples, and knockdown of ETS1 in NA-ATLL cells resulted in inhibition of cell growth. CCR4, a previously identified oncogenic factor in ATLL, was found to be a direct ETS1 transcriptional target in NA-ATLL. As such, ETS1 provides an alternate mechanism to enhance CCR4 expression/activity in NA-ATLL, even in the absence of activating CCR4 mutations (CCR4 mutations were identified in 4 of 9 NA-ATLL cases). Taken together, this study identifies ETS1 as a novel dominant oncogenic transcriptional regulator in NA-ATLL.

A Cellular Stress Response Induced by the CRISPR-dCas9 Activation System Is Not Heritable Through Cell Divisions.

The CRISPR-Cas9 system can be modified to perform "epigenetic editing" by utilizing the catalytically inactive (dead) Cas9 (dCas9) to recruit regulatory proteins to specific genomic locations. In prior studies, epigenetic editing with multimers of the transactivator VP16 and guide RNAs (gRNAs) was found to cause adverse cellular responses. These side effects may confound studies inducing new cellular properties, especially if the cellular responses are maintained through cell divisions-an epigenetic regulatory property. Here, we show how distinct components of this CRISPR-dCas9 activation system, particularly dCas9 with untargeted gRNAs, upregulate genes associated with transcriptional stress, defense response, and regulation of cell death. Our results highlight a previously undetected acute stress response to CRISPR-dCas9 components in human cells, which is transient and not maintained through cell divisions.

Detecting, quantifying, and discriminating the mechanism of mosaic chromosomal aneuploidies using MAD-seq.

Current approaches to detect and characterize mosaic chromosomal aneuploidy are limited by sensitivity, efficiency, cost, or the need to culture cells. We describe the mosaic aneuploidy detection by massively parallel sequencing (MAD-seq) capture assay and the MADSEQ analytical approach that allow low (<10%) levels of mosaicism for chromosomal aneuploidy or regional loss of heterozygosity to be detected, assigned to a meiotic or mitotic origin, and quantified as a proportion of the cells in the sample. We show results from a multi-ethnic MAD-seq (meMAD-seq) capture design that works equally well in populations of diverse racial and ethnic origins and how the MADSEQ analytical approach can be applied to exome or whole-genome sequencing data, revealing previously unrecognized aneuploidy or copy number neutral loss of heterozygosity in samples studied by the 1000 Genomes Project, cell lines from public repositories, and one of the Illumina Platinum Genomes samples. We have made the meMAD-seq capture design and MADSEQ analytical software open for unrestricted use, with the goal that they can be applied in clinical samples to allow new insights into the unrecognized prevalence of mosaic chromosomal aneuploidy in humans and its phenotypic associations.

The Complete Genome Sequence of the Emerging Pathogen Mycobacterium haemophilum Explains Its Unique Culture Requirements.

UNLABELLED: Mycobacterium haemophilum is an emerging pathogen associated with a variety of clinical syndromes, most commonly skin infections in immunocompromised individuals. M. haemophilum exhibits a unique requirement for iron supplementation to support its growth in culture, but the basis for this property and how it may shape pathogenesis is unclear. Using a combination of Illumina, PacBio, and Sanger sequencing, the complete genome sequence of M. haemophilum was determined. Guided by this sequence, experiments were performed to define the basis for the unique growth requirements of M. haemophilum. We found that M. haemophilum, unlike many other mycobacteria, is unable to synthesize iron-binding siderophores known as mycobactins or to utilize ferri-mycobactins to support growth. These differences correlate with the absence of genes associated with mycobactin synthesis, secretion, and uptake. In agreement with the ability of heme to promote growth, we identified genes encoding heme uptake machinery. Consistent with its propensity to infect the skin, we show at the whole-genome level the genetic closeness of M. haemophilum with Mycobacterium leprae, an organism which cannot be cultivated in vitro, and we identify genes uniquely shared by these organisms. Finally, we identify means to express foreign genes in M. haemophilum. These data explain the unique culture requirements for this important pathogen, provide a foundation upon which the genome sequence can be exploited to improve diagnostics and therapeutics, and suggest use of M. haemophilum as a tool to elucidate functions of genes shared with M. leprae. IMPORTANCE: Mycobacterium haemophilum is an emerging pathogen with an unknown natural reservoir that exhibits unique requirements for iron supplementation to grow in vitro. Understanding the basis for this iron requirement is important because it is fundamental to isolation of the organism from clinical samples and environmental sources. Defining the molecular basis for M. haemophilium's growth requirements will also shed new light on mycobacterial strategies to acquire iron and can be exploited to define how differences in such strategies influence pathogenesis. Here, through a combination of sequencing and experimental approaches, we explain the basis for the iron requirement. We further demonstrate the genetic closeness of M. haemophilum and Mycobacterium leprae, the causative agent of leprosy which cannot be cultured in vitro, and we demonstrate methods to genetically manipulate M. haemophilum. These findings pave the way for the use of M. haemophilum as a model to elucidate functions of genes shared with M. leprae.

RNASeq in C. elegans Following Manganese Exposure.

Manganese is a metal that is required for optimal biological functioning of organisms. Absorption, cellular import and export, and excretion of manganese are all tightly regulated. While some genes involved in regulation, such as DMT-1 and ferroportin, are known, it is presumed that many more are involved and as yet unknown. Excessive exposure to manganese, usually in industrial settings such as mining or welding, can lead to neurotoxicity and a condition known as manganism that closely resembles Parkinson's disease. Elucidating transcriptional changes following manganese exposure could lead to the development of biomarkers for exposure. This unit presents a protocol for RNA sequencing in the worm Caenorhabditis elegans to assay for transcriptional changes following exposure to manganese. This protocol is adaptable to any environmental exposure in C. elegans. The protocol results in counts of gene transcripts in control versus exposed conditions and a ranked list of differentially expressed genes for further study.

Transcriptional and metabolic signatures of Arabidopsis responses to chewing damage by an insect herbivore and bacterial infection and the consequences of their interaction.

Plants use multiple interacting signaling systems to identify and respond to biotic stresses. Although it is often assumed that there is specificity in signaling responses to specific pests, this is rarely examined outside of the gene-for-gene relationships of plant-pathogen interactions. In this study, we first compared early events in gene expression and later events in metabolite profiles of Arabidopsis thaliana following attack by either the caterpillar Spodoptera exigua or avirulent (DC3000 avrRpm1) Pseudomonas syringae pv. tomato at three time points. Transcriptional responses of the plant to caterpillar feeding were rapid, occurring within 1 h of feeding, and then decreased at 6 and 24 h. In contrast, plant response to the pathogen was undetectable at 1 h but grew larger and more significant at 6 and 24 h. There was a surprisingly large amount of overlap in jasmonate and salicylate signaling in responses to the insect and pathogen, including levels of gene expression and individual hormones. The caterpillar and pathogen treatments induced different patterns of expression of glucosinolate biosynthesis genes and levels of glucosinolates. This suggests that when specific responses develop, their regulation is complex and best understood by characterizing expression of many genes and metabolites. We then examined the effect of feeding by the caterpillar Spodoptera exigua on Arabidopsis susceptibility to virulent (DC3000) and avirulent (DC3000 avrRpm1) P. syringae pv. tomato, and found that caterpillar feeding enhanced Arabidopsis resistance to the avirulent pathogen and lowered resistance to the virulent strain. We conclude that efforts to improve plant resistance to bacterial pathogens are likely to influence resistance to insects and vice versa. Studies explicitly comparing plant responses to multiple stresses, including the role of elicitors at early time points, are critical to understanding how plants organize responses in natural settings.

Mosaic epigenetic dysregulation of ectodermal cells in autism spectrum disorder.

DNA mutational events are increasingly being identified in autism spectrum disorder (ASD), but the potential additional role of dysregulation of the epigenome in the pathogenesis of the condition remains unclear. The epigenome is of interest as a possible mediator of environmental effects during development, encoding a cellular memory reflected by altered function of progeny cells. Advanced maternal age (AMA) is associated with an increased risk of having a child with ASD for reasons that are not understood. To explore whether AMA involves covert aneuploidy or epigenetic dysregulation leading to ASD in the offspring, we tested a homogeneous ectodermal cell type from 47 individuals with ASD compared with 48 typically developing (TD) controls born to mothers of ≥35 years, using a quantitative genome-wide DNA methylation assay. We show that DNA methylation patterns are dysregulated in ectodermal cells in these individuals, having accounted for confounding effects due to subject age, sex and ancestral haplotype. We did not find mosaic aneuploidy or copy number variability to occur at differentially-methylated regions in these subjects. Of note, the loci with distinctive DNA methylation were found at genes expressed in the brain and encoding protein products significantly enriched for interactions with those produced by known ASD-causing genes, representing a perturbation by epigenomic dysregulation of the same networks compromised by DNA mutational mechanisms. The results indicate the presence of a mosaic subpopulation of epigenetically-dysregulated, ectodermally-derived cells in subjects with ASD. The epigenetic dysregulation observed in these ASD subjects born to older mothers may be associated with aging parental gametes, environmental influences during embryogenesis or could be the consequence of mutations of the chromatin regulatory genes increasingly implicated in ASD. The results indicate that epigenetic dysregulatory mechanisms may complement and interact with DNA mutations in the pathogenesis of the disorder.

Dampened activity of E2F1-DP and Myb-MuvB transcription factors in Drosophila endocycling cells.

The endocycle is a variant cell cycle comprised of alternating gap (G) and DNA synthesis (S) phases (endoreplication) without mitosis (M), which results in DNA polyploidy and large cell size. Endocycles occur widely in nature, but much remains to be learned about the regulation of this modified cell cycle. Here, we compared gene expression profiles of mitotic cycling larval brain and disc cells with the endocycling cells of fat body and salivary gland of the Drosophila larva. The results indicated that many genes that are positively regulated by the heterodimeric E2F1-DP or Myb-MuvB complex transcription factors are expressed at lower levels in endocycling cells. Many of these target genes have functions in M phase, suggesting that dampened E2F1 and Myb activity promote endocycles. Many other E2F1 target genes that are required for DNA replication were also repressed in endocycling cells, an unexpected result given that these cells must duplicate up to thousands of genome copies during each S phase. For some EF2-regulated genes, the lower level of mRNA in endocycling cells resulted in lower protein concentration, whereas for other genes it did not, suggesting a contribution of post-transcriptional regulation. Both knockdown and overexpression of E2F1-DP and Myb-MuvB impaired endocycles, indicating that transcriptional activation and repression must be balanced. Our data suggest that dampened transcriptional activation by E2F1-DP and Myb-MuvB is important to repress mitosis and coordinate the endocycle transcriptional and protein stability oscillators.

Transformation of oats and its application to improving osmotic stress tolerance.

Oat (Avena sativa L.), a worldwide temperate cereal crop, is deficient in tolerance to osmotic stress due to drought and/or salinity. To genetically transform the available commercial oat cultivars, a genotype-independent and efficient regeneration system from shoot apical meristems was developed using four oat cultivars: Prairie, Porter, Ogle, and Pacer. All these oat cultivars generated a genotype-independent in vitro differentiated multiple shoots from shoot apical meristems at a high frequency. Using this system, three oat cultivars were genetically co-transformed with pBY520 (containing hva1 and bar) and pAct1-D (containing gus) using biolistic trade mark bombardment. Transgenic plants were selected and regenerated using herbicide resistance and GUS as a marker. Molecular and biochemical analyses of putative transgenic plants confirmed the co-integration of hva1 and bar genes with a frequency of 100%, and 61.6% of the transgenic plants carried all three genes (hva1, bar and gus). Further analyses of R0, R1, and R2 progenies confirmed stable integration, expression, and Mendalian inheritance for all transgenes. Histochemical analysis of GUS protein in transgenic plants showed a high level of GUS expression in vascular tissues and in the pollen grains of mature flowers. Immunochemical analysis of transgenic plants indicated a constitutive expression of hva1 at all developmental stages. However, the level of HVA1 was higher during the early seedling stages. The characteristic of HVA1 expression for osmotic tolerance in transgenic oat progeny was analyzed in vitro as well as in vivo. Transgenic plants exhibited significantly (P<0.05) increased tolerance to stress conditions than non-transgenic control plants. The symptoms of wilting or death of leaves as observed in 80% of non-transgenic plants due to osmotic stress was delayed and detected only in less than 10% of trans-genic plants. These observations confirmed the characteristic of HVA1 protein as providing or enhancing the osmotic tolerance in transgenic plants against salinity and possible water-deficiency stress conditions.

The Einstein Center for Epigenomics: studying the role of epigenomic dysregulation in human disease.

There is increasing interest in the role of epigenetic and transcriptional dysregulation in the pathogenesis of a range of human diseases, not just in the best-studied example of cancer. It is, however, quite difficult for an individual investigator to perform these studies, as they involve genome-wide molecular assays combined with sophisticated computational analytical approaches of very large datasets that may be generated from various resources and technologies. In 2008, the Albert Einstein College of Medicine in New York, USA established a Center for Epigenomics to facilitate the research programs of its investigators, providing shared resources for genome-wide assays and for data analysis. As a result, several avenues of research are now expanding, with cancer epigenomics being complemented by studies of the epigenomics of infectious disease and a neuroepigenomics program.

Endocycling cells do not apoptose in response to DNA rereplication genotoxic stress.

Initiation of DNA replication at origins more than once per cell cycle results in rereplication and has been implicated in cancer. Here we use Drosophila to examine the checkpoint responses to rereplication in a developmental context. We find that increased Double-parked (Dup), the Drosophila ortholog of Cdt1, results in rereplication and DNA damage. In most cells, this rereplication triggers caspase activation and apoptotic cell death mediated by both p53-dependent and -independent pathways. Elevated Dup also caused DNA damage in endocycling cells, which switch to a G/S cycle during normal development, indicating that rereplication and the endocycling DNA reduplication program are distinct processes. Unexpectedly, however, endocycling cells do not apoptose regardless of tissue type. Our combined evidence suggests that endocycling apoptosis is repressed in part because proapoptotic gene promoters are silenced. Normal endocycling cells had DNA lesions near heterochromatin, which increased after rereplication, explaining why endocycling cells must constantly repress the genotoxic apoptotic response. Our results reveal a novel regulation of apoptosis in development and new insights into the little-understood endocycle. Similar mechanisms may operate during vertebrate development, with implications for cancer predisposition in certain tissues.

Arabidopsis GH3-LIKE DEFENSE GENE 1 is required for accumulation of salicylic acid, activation of defense responses and resistance to Pseudomonas syringae.

In Arabidopsis, the GH3-like gene family consists of 19 members, several of which have been shown to adenylate the plant hormones jasmonic acid, indole acetic acid and salicylic acid (SA). In some cases, this adenylation has been shown to catalyze hormone conjugation to amino acids. Here we report molecular characterization of the GH3-LIKE DEFENSE GENE 1 (GDG1), a member of the GH3-like gene family, and show that GDG1 is an important component of SA-mediated defense against the bacterial pathogen Pseudomonas syringae. Expression of GDG1 is induced earlier and to a higher level in response to avirulent pathogens compared to virulent pathogens. gdg1 null mutants are compromised in several pathogen defense responses, including activation of defense genes and resistance against virulent and avirulent bacterial pathogens. Accumulation of free and glucoside-conjugated SA (SAG) in response to pathogen infection is compromised in gdg1 mutants. All defense-related phenotypes of gdg1 can be rescued by external application of SA, suggesting that gdg1 mutants are defective in the SA-mediated defense pathway(s) and that GDG1 functions upstream of SA. Our results suggest that GDG1 contributes to both basal and resistance gene-mediated inducible defenses against P. syringae (and possibly other pathogens) by playing a critical role in regulating the levels of pathogen-inducible SA. GDG1 is allelic to the PBS3 (avrPphB susceptible) gene.

Overexpression of CRK13, an Arabidopsis cysteine-rich receptor-like kinase, results in enhanced resistance to Pseudomonas syringae.

Protein kinases play important roles in relaying information from perception of a signal to the effector genes in all organisms. Cysteine-rich receptor-like kinases (CRKs) constitute a sub-family of plant receptor-like kinases (RLKs) with more than 40 members that contain the novel C-X8-C-X2-C motif (DUF26) in the extracellular domains. Here we report molecular characterization of one member of this gene family, CRK13. Expression of this gene is induced more quickly and strongly in response to the avirulent compared with the virulent strains of Pseudomonas syringae, and peaks within 4 h after pathogen infection. In response to dexamethasone (DEX) treatment, plants expressing the CRK13 gene from a DEX-inducible promoter exhibited all tested features of pathogen defense activation, including rapid tissue collapse, accumulation of high levels of several defense-related gene transcripts including PR1, PR5 and ICS1, and accumulation of salicylic acid (SA). In addition, these plants suppressed growth of virulent pathogens by about 20-fold compared with the wild-type Col-0. CRK13-conferred pathogen resistance is salicylic acid-dependent. Gene expression analysis using custom cDNA microarrays revealed a remarkable overlap between the expression profiles of the plants overexpressing CRK13 and the plants treated with Pst DC3000 (avrRpm1). Our studies suggest that upregulation of CRK13 leads to hypersensitive response-associated cell death, and induces defense against pathogens by causing increased accumulation of salicylic acid.